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ABSTRACT 

We present a pseudo-Newtonian potential for accretion disk modeling around the rotating black holes. 
This potential can describe the general relativistic effects on accretion disk. As the inclusion of rotation 
in a proper way is very important at an inner edge of disk the potential is derived from the Kerr metric. 
This potential can reproduce all the essential properties of general relativity within 10% error even for 
rapidly rotating black holes. 
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1. INTRODUCTION 

Most of the theoretical studies of general relativity in 
astronomy are approached by a Newtonian or pseudo- 
Newtonian method. To avoid the complexity of full general 
relativistic equations it is simpler to use non-relativistic 
equations but with the inclusion of corresponding (pseudo) 
potential which can reproduce some relativistic effects ac- 
cording to the geometry of space-time. Using this poten- 
tial one can get the approximate solutions of the hydro- 
dynamical equations. Shakura & Sunyaev (1973) initiated 
the modeling of accretion disk around black holes using 
Newtonian gravitational potential as 



(1) 



where, x = r/rg, r is the radial coordinate of the disk and 
Tg = GM/c^. However, this potential can not reproduce 
the properties of the inner region of a disk where relativis- 
tic effects become important. Later on, Paczyhski & Wiita 
(1980) proposed a pseudo-Newtonian potential which can 
reproduce approximately the properties of the inner disk 
at close to the equatorial plane around non-rotating black 
holes without using relativistic fluid equations as 



(x-2)- 



(2) 



The beauty of their potential is that, it can give the right 
positions of the marginally stable (xs) and marginally 
bound (xb) orbits in a Schwarzschild metric. It also re- 
produces the total mechanical energy per unit mass at the 
last stable circular orbit (Eg) and the total energy dissipa- 
tion at a given radius (77) in good agreement with that of 
Schwarzschild geometry (Artemova et al. 1996). The error 
in both cases is less than 10%. Nowak & Wagoner (1991) 
proposed another potential for accretion disk around non- 
rotating black holes as 
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which can reproduce the values of Xs and angular velocity 
(fi) at that radius in Schwarzschild geometry. So far this 



choice of the potential gives the best approximate radial 
epicyclic frequency. Artemova et al. (1996) proposed two 
correct potentials for describing the accretion disk around 
rotating black holes. The form of one of their potentials 
is given as 



dV3 
dx 
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K^-P{x-xi)P' 
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xi is the black hole horizon and /3 is a constant for a par- 
ticular specific angular momentum of the black hole, a (for 
exact expression see Artemova et al. (1996)). They showed 
that their potentials can reproduce the value of Xs exactly 
as that for Kerr geometry and reproduce the values of r] at 
different radii in good agreement with that of general rela- 
tivistic results. After that, several authors (e.g. Artemova 
et al. (1996), Lovas (1998), Semerak & Karas (1999)) have 
analyzed the efficiency of different pseudo-potentials pre- 
scribed for accretion disks around black holes. In the con- 
text of accretion-disk-corona, which is infalling towards 
the black hole, Miwa et al. (1998) chose the pseudo- 
potential V3 (Eqn. (^) and discussed about radiation flux, 
velocity of infalling corona upto very close to the black 
hole. 

However, while Eqn. (H) is the analytical form for force, 
other pseudo-potentials Tor non-rotating black holes like 
Eqns. (^) and ^ have a simple analytical form of poten- 
tial. For the study of parameter space, e.g. sonic point 
analysis etc. it is useful to have an analytical expression 
for the potential which should asymptotically vary as—l/x 
and reduce to Paczyhsky- Wiita form (Eqn. (0)) for zero 
rotation. Apart from that, equation (15) of Artemova et 
al. (1996) is only valid for co-rotating (positive values of 
Kerr parameter, a) disks. If one takes negative values of 
a to represent counter-rotation and using that equation 
(15) calculates Eg and Xf,, the error may be upto 50% 
and 500% respectively. Similarly, for negative a it can 
not give the correct value of Xg- The general expression 
of equation (15) of Artemova et al. (1996) should read 
Hn 3 -h ^2 T [(3 - ^i)(3 + Zi + 2Z2)]i/2 (Bardeen 1973, 
Novikov & Thorne 1973), where the upper and lower signs 
are for co-rotation and counter-rotation respectively. 

As the inner region of an accretion disk is very influ- 
enced by the rotation of black hole, rotation should be in- 
corporated in theoretical studies. There are observational 
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indications that black holes could be rotating rapidly and 
thus to study the inner properties of an accretion disk ro- 
tation should be incorporated correctly. Iwasawa et al. 
(1996) have argued that the variable iron K emission line 
in MCG-6-30-15 arises from the inner part of an accretion 
disk and it is strongly related to the spin of the black hole. 
It has been argued from other observational point of view 
(Karas & Kraus 1996, Iwasawa et al. 1996) that central 
black holes in galactic nuclei is likely to be rapidly rotating. 
Also for the observation of gravitomagnetic precession, the 
inner edge of accretion disk is responsible (Markovic & 
Lamb 1998, Stella & Vietri 1998). The temporal prop- 
erties of the system are expected to depend on the inner 
edge of disk which in turn depends on the rotation of the 
black hole. The predictions of the disk properties will be 
incorrect if the pseudo-Newtonian modeling does not take 
into account spin of the black hole. 

The aim of this paper is to present a pseudo-Newtonian 
potential which can reproduce exactly or in good agree- 
ment all the (inner) accretion disk properties close to the 
equatorial plane in Kerr geometry. The potential should 
reproduce those features of a rotating black hole geometry 
which have been reproduced by Paczyiiski & Wiita (1980) 
potential (Eqn. ||) for a non-rotating black hole. Thus, we 
will establish our potential in a same spirit as Paczyhski 
and Wiita did for a non-rotating black hole. All other 
forms of the potential have been introduced without clear 
relation to the space-time metric. Here we will formulate 
our pseudo-potential from the Kerr metric. As the metric 
is involved directly to our calculation many of the features 
of Kerr geometry are inherent in our potential by design. 
In the next section we present the basic equations and de- 
rive the pseudo-potential. In §3, we compare a few results 
of the Kerr geometry with that of the potential and in §4, 
we make our conclusions. 

2. BASIC EQUATIONS AND PSEUDO-POTENTIAL 

The Lagrangian density for a particle in the Kerr space- 
time in Boyer-Lindquist coordinate at the equatorial plane 
{9 = 7r/2) can be written as 



2C = -[l- 



2GM\ AGMa ■ ■ 



t(j) + — 



a2 2GA/a2 



(5) 



where over-dots denote the derivative with respect to the 
proper-time r and IS. — r"^ + c? jt? — 2GMr/c^. 



The geodesic equations of motion are 



E — constant = 



2GM\ . 2GMa 
1 —]t + —, c 



c^r J 
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2GMa. 

A = constant = 1 - 



2GMa' 



(7) 



For the particle with non-zero rest mass g^i^p^p'^ — ~m? 
(where is the momentum of the particles and g^^ is the 
metric). Replacing the solution for i and from (^ and 



(|^) into (||) gives a differential equation for 
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Here, 5* can be identified as an effective potential for the 
radial geodesic motion. The conditions for circular orbits 
are 



* = 0, 



dr 



= 0. 



Solving for E and A from (bf) we get 



E 



- 2GMr/c^ + a^GMr/d^ 



and 



m r{r^ - 3GMr/c^ + 2a^ GMr j c^ fl'^ 



A _ ^GMrlc^{r^ - 2a^GMr/c^ + a^/c^) 
m ~ r{r^ - 3GMr/c^ + 2a^GMr j c'^yl'^ 



(9) 



(10) 



(11) 



Equations (10) and (11) have been derived by Bardeen 
(1973). Now as standard practice, we can define the Ke- 



plerian angular momentum distribution \k 



There- 



fore, corresponding centrifugal force in Kerr geometry can 
be written as 



X2 



2\2 



(x^ - 2aV^+ a^) 
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Thus from above, Fx can be identified as the gravitational 
force of black hole at the Keplerian orbit. The above 
expression reduces to Paczynski- Wiita form for a = 0. 
Thus we propose Eqn. dl2) is the most general form 
of the gravitational force corresponding to the pseudo- 
potential in accretion disk around black holes. The gen- 
eral form of the corresponding pseudo-potential (which is 
Vx = V4 = J Fxdx) is algebraically complicated, but sim- 
plifies for any given value of a. In Appendix the general 
form of the potential and its reduced form for a few par- 
ticular Kerr parameters are given. 

3. COMPARISON OF THE RESULTS FOR KERR GEOMETRY 
AND PSEUDO-POTENTIAL 

To establish the validity of this potential, we raise the 
following questions. (1) Does this potential {V4) repro- 
duce the values of Xh and Xg as same as Kerr geometry? 
(2) Does it give the correct value of Eg as around of Kerr 
black hole? (3) How does the corresponding dissipation 
energy distribution ri(x) in the accretion disk by this po- 
tential match with that of pure general relativistic result? 

Apart from that, one can ask how simple the form of this 
potential so that it is applicable for other studies where an 
analytical form is required? Below, we are discussing all 
the questions one by one. If we can show, our potential 
tackle all the above issues fairly well, we can conclude that 
this is one of the best potential for an accretion disk around 
rotating back holes as well as non-rotating ones. 

At the marginally bound orbit mechanical energy E re- 
duces to zero and we get 



(13) 
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which can be used to calculate for t4- For the stability 
of an orbit d\/dx > 0, which for the specific potential V4 
is 

- 3a* + 14a^^/^ + {x - 6)x^ + 6ax^^^{x + 2) 

~2a^x{x + ll)>0. (14) 

The solution of Eqn. (ij) with 'equals to' sign gives the lo- 
cation of the last stable circular orbit (xg) for V4. For any 
a, the Xg computed from the above equation dl4) matches 
exactly with the radius of last stable circular orbit in Kerr 
geometry. We are not reporting the Xs for various a values 
as it is available in standard literatures. In Table- 1 and 2 
we list Xb and Eg for the potential V4 and Kerr geometry 
for various values of a. 

Table- 1 
Values of Xb 



a 





0.1 


0.3 


0.5 


0.7 


0.998 


Vi 


4.0 


3.788 


3.347 


2.870 


2.333 


1.037 


Kerr 


4.0 


3.797 


3.373 


2.914 


2.395 


1.091 


a 





-0.1 


-0.3 


-0.5 


-0.7 


-0.998 


Va 


4.0 


4.206 


4.606 


4.993 


5.368 


5.911 


Kerr 


4.0 


4.198 


4.580 


4.949 


5.308 


5.825 



From Table-1, it is clear that for all values of o, V4 can 
reproduce the value of xjj in very good agreement with 
general relativistic results. The maximum error in xi, is 
^ 5%. Table-2 indicates that V4 produces Eg in a fairly 
good agreement with Kerr geometry with a maximum pos- 
sible error ~ 10%. Thus, the potential V4 will product a 
slightly larger luminosity than the general relativistic one 
in the accretion disk for a particular accretion rate. Note 
that, for counter-rotating black holes the errors are less 
than those of a co-rotating one. 

Next we will compare the total energy dissipation 77 
in the accretion disk for this potential with that of gen- 
eral relativistic result. We choose a simple a-disk model 
to compute pseudo-Newtonian 77(2;) (Bjornsson & Svens- 
son 1991, Frank et al. 1992, Artemova 1996), where for 
different choices of pseudo-potential, 51 can be different 
which could have been reflected to the final profile of ri{x) 
(Bjornsson & Svensson 1991, 1992). Following Novikov 
& Thorne (1973), Page & Thorne (1974) and Bjornsson 
(1995) we calculate the corresponding general relativistic 
77 profile. In Fig. 1, we show •q/m (m is the accretion rate) 
as a function of x, for different values of Kerr parameter 
a. Here also our pseudo-Newtonian results agree within 
10% of general relativistic values. For moderate rotation 



almost there is no error, while for very rapidly rotating 
black holes the deviation increases but still within 10%. 

Apart from that, the analytical expression for the force 
(Eqn . (|l2[) ) as well as potential at given values of a (Eqns. 
(A3) and (A4)) are relatively simple which will make it 
easy to implement in other applications like detailed fluid 
dynamical studies, analysis of the parameter space in disk 
etc. Thus this potential satisfies all the criteria for a good 
pseudo-potential which can describe an accretion disk us- 
ing non- relativistic equations. 

4. CONCLUSIONS 

We have prescribed a general pseudo-potential for the 
modeling of accretion disks around rotating black holes. 
For the non-rotating case, it reduces to Paczyhski-Wiita 
potential. Unlike previous works, this potential is derived 
from the metric (Kerr geometry) at the equatorial plane. 
Naturally it exhibits better accuracy as it is derived from 
the metric itself. Following the same procedure, using 
Schwarzschild metric one can derive Paczyhski-Wiita po- 
tential. The detailed calculations of various geodesic equa- 
tion are available in standard literatures (e.g. Shapiro 
& Teukolsky 1983). Our potential is valid for both co- 
rotating and counter-rotating black holes which is not nec- 
essarily the case for earlier potentials (Artemova et al. 
1996). In fact for our pseudo-Newtonian potential, the 
counter-rotating results agree better with the general rel- 
ativistic ones, presumably because of the larger values of 
Xb and Xg with respect to that of co-rotating cases. But, 
still for our potential the possible error is 10% at most for 
any rotation. Thus the inference of various observational 
aspects using our potential may have better accuracy. If 
the description of disk property is acceptable within 10% 
accuracy, our potential should be recommended. It should 
be mentioned that, although this pseudo-Newtonian po- 
tential is applicable close to the equatorial plane, it may 
not be a good approximation to follow light rays, orbits 
far from the equator. As because at the very beginning 
of our calculation we have chosen = n/2 this further 
constraint arises. Such a pseudo-Newtonian potential for 
generalized 6 has not been proposed before and hence it 
might be useful to derive that one following the method 
prescribed in this work. 

Next, one can apply this potential for various fluid dy- 
namical problems. It is shown here that the maximum er- 
ror in calculation of r] is 10% even for the rapidly rotating 
black holes. One should study: how does the rotation of 
black hole affect the fluid properties in an accretion disk? 
How does it affect the parameter region of disk? In a next 
work, we expect to explore all these issues. 



APPENDIX 

APPENDIX: ANALYTICAL EXPRESSIONS OF PSEUDO-POTENTIAL 

As we have an analytical form of the gravitational force ([l^), we can calculate the corresponding potential as Vx 
J Fxdx. Thus the most general expression for the pseudo-potential is 

4a 2{%a^x-lQax+l&^/x-lia^^/x + &a?^ -8a) 
*^""2^+ (27a2-32)(x3/2-2V^ + a) ^^^""^ 
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Table-2 
Values of Eg 



(1 







0.1 


().:-! 


0.5 


0.7 


0.998 


Vi 




0.0625 


-0.0663 


-0.0761 


-0.0904 


-0.1149 


-0.3533 


K(^rr 




0.0571 


-0.0606 


-0.0693 


-0.0821 


-0.1036 


-0.3209 









-0.1 


-0.3 


-0.5 


-0.7 


-0.998 


Vi 




0.0625 


-0.0591 


-0.0536 


-0.0491 


-0.0454 


-0.0409 


Kerr 




0.0571 


-0.0542 


-0.0492 


-0.0451 


-0.0418 


-0.0378 




Fig. 1. — Energy dissipation per unit accretion rate (ri/m) as a function of radial coordinate {x) for various values of specific angular 
momentum of the black hole (o), which are indicated at each set. For each set, solid and dotted curves show results of general relativity and 
our potential respectively. 
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27a2 - 32 
74a?/ - 107a 



E 



y=Xi.X2,X3 

2 



-logiV^ - y){54a^y^ ~ 6iy 



where, 



Xl 



24/3 

p 



128)] , 



21/33' 



(Al) 



X2 



21/36 



and p= (v/729a2 - 864- 27a) 



1/3 



(l + iV3). 



(A2) 



Here, denotes the complex conjugate. Though Eqn. (Al) looks very complicated if we specify the particular values 
of a it reduces to rather simpler expression. For example, if we choose a = 0, it reduces to Paczynski-Wiita potential 
(— j^). For some other values of a the analytical forms of the potential reduce as 



ya=±l 

where, A = 2.15542, B 



2x2 
1.61803, C 



1 4 
± 



h loq - — = -77 



(a;3/2-2Vi±l) 
12.1554, D = 0.618034 and 



2log{x), 



(A3) 



a=±0.5 



1 _^ 2 ^ 
Ba;2 -Vx 



T7.75x + 25.5ViT6.5 



log 



2x3/2 -4V^± 1 



2log{x), 



,(V^T/)^(V^Ti^)'' 
1.52569, iJ = 5.93863, / = 1.26704, J = 50.2075, K = 0.258652. As a 



(A4) 
^ 1, 



where, E = 0.0792079, = 5.64616, G 
iJ — > 0. One can easily check that both (A3) and (A4) asymptotically vary as — 1/x. The expressions (A3) and (A4) 
could have been reduced to a more simpler form if we approximate the values of decimal number in constants. But in 
this manner accuracy of the solution would be reduced. 
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